BKT phase transition in a 2d system with long range dipole-dipole interaction. 
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We consider phase transitions in 2d XY-like systems with long range dipole-dipole interactions and 
demonstrate that BKT-type phase transition always occurs separating the ordered (ferroelectric) 
and the disordered (paraelectric) phases. The low-temperature phase corresponds to a thermal 
state with bound vortex- ant i vortex pairs characterized by linear attraction at large distances. We 
show that bound vortex pairs polarize and screen the vortex-antivortex interaction leaving only the 
logarithmic attraction at sufficiently large separations between the vortices. At higher temperatures 
the pairs dissociate and the phase transition similar to BKT occurs, though at a larger temperature 
than in a system without the dipole-dipole interaction. 



Phase transitions and long range order in 2d systems 
are intimately related with dynamics of topological exci- 
tations (vortices) with logarithmic interactions: the tran- 
sition occurs when vortex-antivortex pairs dissociate at 
a certain finite temperature[ll, [3, 0]. This conclusion is 
well estabHshed for a free system and does not quali- 
tatively change in a presence of a short range interac- 
tion. The ground state of the system is not truly ordered 
(there can be no true long range order) but superfluid 
at temperatures below the phase transition point. The 
transition has been observed in numerous experimental 
systems, such as, most recently, weakly interacting ultra- 
cold gases [3|. 

There are situations though when the microscopic in- 
teractions are long range, as in a case of a system with 
appreciable dipole-dipole forces. The nature of the phase 
transition under the circumstances has a long history 
[l,0,0,[l,S[l3,[lH- Dipole-dipole interaction are shown 
to lead to formation of a true long range order at small 
temperatures, so that the ground state of a 2d system of 
magnetic dipoles is a ferromagnet The spontaneous 
polarization of the ground state has been also pointed out 
in [fl] . Both works do not clear out the nature of the phase 
transition, which has been studied first in where the 
correction to BKT temperature was obtained in the limit 
when the dipole-dipole interaction is much smaller than 
the exchange term, and hence the correction to the tran- 
sition temperature is small. The author made an impor- 
tant observation though: the interaction between the vor- 
tices was found to be linear at sufficiently large distances 
(see for a recent discussion) . Ferromagnetic ordering 
was studied in various model dipole systems theoretically 
in 0, H, [§| and experimentally in [llj . Arguably, dipole- 
dipole forces between water molecules on hydrophobic 
surfaces lead to macroscopic ordering of the molecular 
dipoles such as creation of macroscopic hydrogen-bond 
networks in biological systems pJi]. 

In spite of vast volume of the research, the nature of the 
disordering phase transition appears to be somewhat con- 



troversial. There are claims that the long range dipole- 
dipole forces change the physics of the phase transition 
entirely transforming BKT vortex pairs dissociation tran- 
sition in a free system to a confinement transition similar 
to that in quark-gluon plasma Note, that transition 
temperature is predicted to be four times higher than 
BKT temperature. Below we perform a systematic study 
of dipole-dipole interaction infiuence in 2d systems at fi- 
nite temperatures. We show that the interaction leads to 
ferromagnetic ordering at low temperatures without con- 
tradiction to Mermin-Peierls theorem and in a full 
accordance to the earlier statements. The phase transi- 
tion itself turns out, as in BKT case, to be associated 
with dissociation of vortices. We calculate the transition 
temperature, Tc, as a function of the interactions param- 
eters. We show that the Hnear interaction between the 
vortices does change the transition temperature, though 
the vortex gas polarization screens the long range linear 
potential and transforms it to a logarithmic interaction. 
Therefore, in a system with dipole-dipole interaction the 
phase transition is essentially BKT, though the transi- 
tion temperature itself has a complicated dependence on 
the model parameters. 

At last we apply the developed model to water-solute 
boundaries using the phenomenological vector model of 
polar Hquids We show that spontaneous polariza- 

tion of molecular dipoles next to hydrophobic boundaries 
may occur and a ferroelectric liquid film may form. The 
vector model naturally describes topological excitations 
on the solute boundaries. Since the molecular dipoles in- 
teract electrostatically at large distances, the dissociation 
of the vortex pairs can be associated with disappearance 
of hydrogen bond networks in the course of order-disorder 
phase transition in the hydration water layer [15, 16, 17J. 

Consider a plane layer of a thickness A composed of 
interacting dipoles. The unit vector S(r) is taken parallel 
to the dipole moment of a molecule residing at a point r 
and is characterized by the orientation angle 6{t): S = 
(cos^^(r), sin0(r)). The Hamiltonian of the interacting 
system consists of the two parts: 
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Gs — Gh + Gdd, (1) 

where the short ranged gradient term gives the energy of 



2 



a free system 



Gh = \m I df {vef 



(2) 



and the long range dipole-dipole interaction is repre- 
sented by the following model term: 



G, 



dd 
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dfdf 



,(V.S) (V'-S') 



where df is the element of the surface F, whereas the con- 
stants M and K characterize the strength of the inter- 
actions. The Hamiltonian has been studied in the limit 
M:^ XK in [i,[i3|. 

Consider the case of an arbitrary relation between pa- 
rameters XK and M. At small temperatures T the dipole 
system is ordered, i.e. all the dipoles point at the same 
direction, for example, along the x-axis {9{r) =0). To 
prove that let us consider the correlation function of 
the dipole orientation fluctuations: K2 (r) = {9 (0) 6 (r)), 
\9\ <C 1. Keeping quadratic terms in 9 only and lineariz- 
ing the Hamiltonian we find that 
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Gs ~ t;M / df {^9Y + -K / dfdf 



ir)9y (r') 



so that the correlation function takes the form 

exp (ikr) 
Mfc2 + 27rKk sin^ a ' 



(3) 



where a is the angle between wave vector k and a;-axes. 
The integral in the r.h.s. converges at small values of 
k, therefore i^2(0) is finite and hence the ordered state 
is thermally stable. In fact the appearance of the linear 
term in the integrand denominator resolves the contradic- 
tion with Mermin-Peierls theorem It may even seem 
that the correlation function integral converges solely due 
to the finite value of the dipole dipole interaction strength 
{K ^ 0). The conclusion is in fact wrong: as it is clear 
from Eq. ([3]), at M ^ the integral diverges at a — > 0. 
It means that both terms in Eq. ([!]) denominator are 
both equally important. Mathematically this analysis is 
in full agreement with a well known feature of classic 
electrostatics, the Earnshaw theorem: a system of classi- 
cal charges interacting with electrostatic forces only can 
have no stable state 

To elucidate the nature of the ordered state let us con- 
sider the correlation function at large distances r. The 
main contribution to the integral in Eq. |(3| comes from 
the two separate regions: |q;| ^ 1 and |a — 7r| ^ 1, so 
that 



K2{v) 



4T , 
-r-F= I dr cos 
Jo ^ 



XT 1 exp 



^/7 



(4) 



with 7 = 2i:K/M. The calculation in Eq.|[4l) simplifies 
in the two following limiting cases: 
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Figure 1: Polarization configuration S{x,y) corresponding to 
a vortex- ant i vortex pair (VAP). 



The correlations at large temperatures are addressed in 
[lo| and turn out to decay exponentially. The crossover 
between the power laws found here for small T and the 
exponential decay at large T suggests a phase transition 
at some intermediate temperature Tc- 

To establish the transition temperature we apply a 
set of familiar arguments (see e.g. j3|| for further 
explanations). Thermodynamics of a 2d system can 
be mapped to thermodynamics of a gas of interacting 
vortex-antivortex pairs (VAPs) [l^. Consider first the 
small temperatures limit :T < Tc, Tc — T <^ Tc. The 
thermal state in this case can be viewed as a gas of bound 
VAPs. Since the considered temperature is still below the 
transition point, VAPs approach the dissociation limit 
and the pairs with large distances r between the vortices 
dominate. At such separations the vortex-antivortex in- 
teraction is described by the linear term and the energy 
of the system can be approximated as: 



F^NEo + \k J2 QrQj (.-rrj) = J E (^«) 



(5) 



where N is the number of the VAPs and qj — ±1 are 
the topological charges associated with the vortices. The 
"rest energy", Eq, is the energy associated with a pair 
of a minimum possible separation r ~ A, where A is the 
size of the vortex core. The latter quantity is is small, 
model dependent and depends on microscopic details of 
the underlying physical system (such as thickness of the 
layer). The (quasi-electric) potential <& (r) is introduced 
in analogy with classic electrostatics and amounts to the 
energy associated with the interaction of a given vortex 
(charge) with all other vortices in the system: 
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(6) 
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where 



>{r)=J2q,6^'^ (r-r,) 



(7) 



is the vortex charge density. Note, that according to its 
definition ([6]) , the potential $ (r) satisfies the analogue 
of Poisson equation 

Lr$(r) =ifp(r), 
where the linear operator is defined so that 

LA-\r-r'\) = S<^'Hv-r'). (8) 

Let us follow the electrostatic analogy even further: since 
the energy of a charge (vortex) q placed in the external 
potential $ is J7 = q^, then the force, acting on the 
charge (vortex) is F = qE, where the quasi-electric field 
vector E = — V$ is associated with the potential $ in a 
normal way. 

To find out the energy of a VAP in a self-consistent 
"electric" field of all other pairs, let us define the "dipole 
moment" of a pair according to 



(9) 



where r = r+ — r_. Then, the energy of the pair is 



U = (Z+$(r+) + g_$(r_; 



d E. 



(10) 



Let us pursue the analogy and calculate first the polariz- 
ability ap of a single VAP. The dipole moment of a pair 
in a weak external field E is given by a standard relation 



J dfr exp 



r-E 

T 



K„ 



/rf/exp(^^-fr 



apE, 



(11) 



where ap = 3T/K^ is nothing else but the pair polariz- 
ability. 

At the transition temperature, T = Tc, VAPs begin 
to dissociate. It means that at Tc — T <C lb only a 
small fraction of the pairs are very large and close to 
dissociation. For this reason it is possible to neglect the 
interactions between the largest VAPs and calculate the 
energy of single large pair approaching its dissociation 
limit in a cloud of comparatively small bound VAPs. As 
we have demonstrated above, the bound pairs are polar- 
izable and therefore the field of a charge is screened by 
the polarization of VAPs gas, thus infiuencing the poten- 
tial energy of a large VAP. To find out how exactly, let us 
consider the effect of shielding of a probe point charge Q 
or the charge density pg — QS'^^^ (r) placed somewhere at 
the origin. The complete "electrostatic" potential $ (r) 
is produced both by the probe charge and the polariza- 
tion charges of VAPs. The density pp of the polarization 
charges is given by usual expressions following from Eq. 
([9]): pp = — V • P, where P is polarization of the vortex 



gas: 



np (d) = xE 



(12) 



with X = ctpnp and np being the "dielectric suscepti- 
bility" of the gas of VAPs and the concentration of the 
pairs, correspondingly. 

Combining Eqs. ([8]), (fT2|) and ([6]) we obtain the com- 
plete equation for $ (r) : 



(r) ^K{pQ + pp) = KQS<-^^ (r) + if xA$ (r) 



(13) 



Since the Fourier component of the operator L^. is ik 
k^/2Tr, we find that 



$(r) = QK 



(Pk exp (ikr) 



(14) 



At large distances, r ^ rg — 1/2-kKx, most contribution 
to the integral comes from small values of k where k^ term 
in the denominator is negligible. Therefore the potential 
of the "charge" at large distances is logarithmic: 



(15) 



where Ci ~ tq. In the opposite limit, r <C rp, the poten- 
tial is Hnear: $(r) « —Kr, r <C tq. Below we propose a 
simple expression interpolating between the two results: 



(16) 



27rx V '"0 



so that the energy of a large pair of a size R is given by 

(17) 



1 I R 

EvAp[R) « 2i;o + — log 1 + — 
27rx V ^0 



We note that although the dipole-dipole interactions do 
change the interactions between the vortices at small dis- 
tances, the polarization of the VAPs destroys the Hnear 
attraction at large separations between the vortices and 
hence the phase transition associated with the dissoci- 
ation of the pairs is qualitatively very similar to BKT 
transition in a free system. 

Standard calculation of BKT temperature for a vor- 
tex gas with interaction (fTTl gives the following implicit 
equation for the transition temperature Tc 



12Tmp(Tc) 



(18) 



The density of np at the transition temperature can be 
calculated from the following argument. At T « Tc the 
VAPs just start to dissociate and therefore the size of 
a typical pair is small and the interaction between the 
vortices is still linear and unscreened and 



np(T) ~ 27rexp 



2^0 \ 



T2 



(19) 



Therefore, the Eq. (fT8|) for the transition temperature 
takes form 



En 



(20) 
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Figure 2: Graphical solution of implicit equation l[20|l for the 
BKT-transition temperature Tc- By R and L graphics of the 
r.h.s. and l.h.s. are shown, correspondingly. 

The graphical solution of the equation is demonstrated 
on Fig. [21 As it is clear from the graph, the solution of 
the equation always exists and is unique. 

Let us assume first that the dipole-dipole interactions 
is small, i.e. XK <^ M. Since the energy of a minimal 
vortex- antivortex pair Eq ^ M, the transition temper- 
ature Tc ~ M, which coincides with the standard BKT 
result (Tbkt — 7rM/2). This is precisely the limit con- 
sidered in In the opposite limit, M <^ XK, the 
exponent in the r.h.s. of Eq. (|2Q|) is approximately 1 
and hence Tc '■^ XK 3> Tbkt- The phase transition 
temperature is indeed larger than BKT temperature, in 
agreement with although our analysis suggests that 
Tc can be arbitrary larger than Tbkt- 

Layers of water molecules on hydrophobic surfaces such 
as biological membranes or large bio-molecules are an- 
other interesting example of 2d systems with dipole- 
dipole interactions. The orientations of molecules next 
to a macromolecular surface depend on interplay of long- 
range interactions of the molecular dipoles and the short- 



range hydrogen bonds between the adjacent molecules. 
Since an energy of uncompensated hydrogen bond is large 
compared to the temperature, no water molecules can 
point in the direction of a hydrophobic surface, and thus 
all the water molecules next to the surface have their 
dipole moments parallel to the surface. Macroscopic po- 
larization of the molecules s (r) = (S (r)) vanishes quickly 
in the direction of the liquid bulk. Therefore a water layer 
next to a hydrophobic boundary can be studied with a 
help of a model Hamiltonian Jl]) corresponding to a 2d of 
interacting dipoles (see fl3| and references therein). Ac- 
cording to the findings of the current work the layer of 
water molecules is completely polarized, s « S||, at very 
low temperatures. Microscopic parameters of water are 
such that M ~ XK and therefore the macroscopic polar- 
ization and molecules orientation ordering disappears at 
temperature 

Tc^M^EhsI (21) 

where Eh ~ 2500°^ is the characteristic energy of a 
hydrogen bond. The expression suggests that the transi- 
tion temperature corresponding a completely hydropho- 
bic surface (s|| w 1) is always very large, the water 
molecules are ordered at all realistic temperatures and 
thus there always exists a macroscopic hydrogen bonds 
network on the surface of the body. If a surface is partly 
hydrophillic, i.e. hydrogen bonds donors or acceptors , 
such as charges, sy < 1 and already at S|| — 0.3 the tran- 
sition may occur already at room temperatures. The 
evidence of such transitions is observed in molecular dy- 
namics calculations [H, [H, [l3| ■ 
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ticals. Phenomenological water models are used in 
Quantum everywhere to compute solvation energies and 
protein-drug interactions in aqueous environment. 
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